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We show how to optimize the process of high-harmonic generation (HHG) by gating the inter- 
action using the field gradient of the driving pulse. Since maximized field gradients are efficiently 
generated by self-steepening processes, we first present a generalized theory of optical carrier-wave 
self-steepened (CSS) pulses. This goes beyond existing treatments, which only consider third-order 
nonlinearity, and has the advantage of describing pulses whose wave forms have a range of sym- 
metry properties. Although a fertile field for theoretical work, CSS pulses are difficult to realize 
experimentally because of the deleterious effect of dispersion. We therefore consider synthesizing 
CSS-like profiles using a suitably phased sub-set of the harmonics present in a true CSS wave form. 
Using standard theoretical models of HHG, we show that the presence of gradient-maximized re- 
gions on the wave forms can raise the spectral cut-off and so yield shorter attosecond pulses. We 
study how the quality of the attosecond bursts created by spectral filtering depends on the number 
of harmonics included in the driving pulse. 



I. INTRODUCTION 

In high-order harmonic generation (HHG), the inten- 
sity maxima of a strong driving laser pulse tunnel-ionize 
an atom or molecule. Once an electron is in the con- 
tinuum, its acceleration is controlled by the electric field 
profile E(t). The electron is initially driven outwards, but 
then the field reverses, driving it back for a high-energy 
recollision with the core, approximately three-quarters of 
a cycle later [Q, |) , in which an attosecond burst of ex- 
treme ultraviolet (XUV) radiation is generated. Indeed, 
each half-cycle of the driving pulse initiates this sequence 
of events, so a train of XUV bursts a few hundred attosec- 
onds long is produced. While a train is useful for probing 
ultrafast dynamics, a primary goal of attosecond science 
is the generation of single bursts of XUV radiation, and 
various gating schemes have been employed to achieve 
this. One scheme already in use is polarization gating 
j| [§, ||, ||; in this technique, the sense of circular po- 
larization reverses smoothly through the driving pulse, 
passing only briefly through linear polarization at which 
point an attosecond burst is generated. 

Another technique is to "gradient-gate" the HHG in- 
teraction, which can be done most simply by an ad- 
ditional driving field at double the frequency, suitably 
phased with respect to the fundamental field; this is just 
the case of two-color pumping J?], [§[ [| [To) . In the present 
paper, we push the concept of gradient-gating to its lim- 
its by considering how to achieve the steepest gradients 
possible. In principle, this can be done with carrier- 
wave steepening and shocking, a process first identified 
by Rosen |LlJ in 1965. He showed that for a instantaneous 
third-order (x ) nonlinearity, field gradients could be- 
come infinite (i.e. shocks could form) after a finite dis- 
tance of propagation. As we demonstrate in Section ||, 



the theory can be extended to include arbitrary combina- 
tions of nonlinear terms, which opens up possibilities for 
controlling the symmetry of the pulse at the same time 
as generating a gradient-gate. We then proceed to dis- 
cuss the previously unaddressed second-order x case, 
which turns out to give wave forms particularly useful for 
HHG. 

Unfortunately it is very difficult to generate carrier- 
wave self-steepened (CSS) pulses because dispersion over- 
whelms the steepening process In Section III, we 



therefore discuss the possibility of synthesizing the steep- 
ened profiles by combining a suitably phased sub-set of 
the harmonics present in a true CSS profile; optical wave 
form synthesis of this kind is already exploited in other 
contexts (see, e.g., 0). Such a synthesis is in any 
case the logical extension of both the two-color pump- 
ing scheme and the gradient-gate concept. 

In Section IV, we compare simulations of HHG driven 
by standard pulses and by synthesized CSS-like wave 
forms. In the latter case, we demonstrate the importance 
of the relative phasing of the spectral components by 
comparing the results to those based on phase-flattened 
(i.e. field-maximized) pulses with equivalent spectral 
content. We also show that the efficacy of the synthesized 
profiles for HHG is progressively enhanced as the num- 
ber of harmonics is increased. High-pass filtering then 
leads to the production of attosecond XUV bursts, and 
the presence of even harmonics in x^ CSS-like pulses is 
shown to increase the degree of isolation of the bursts. A 
genetic algorithm is used to show that sawtooth profiles 
are ideal for driving HHG, which is in line with the gen- 
eral conclusions of this study, as summarised in Section 
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The self-steepening of an optical pulse envelope was 
first studied by DeMartini et al. in 1967 (h|], and is a 
well-known phenomenon associated with self-phase mod- 
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ulation. Surprisingly, however, the possibility of shock 
formation on the optical carrier was considered even ear- 
lier in a 1965 paper by Rosen p| , who showed that, 
for a third order (x ) nonlinearity, a field discontinu- 
ity (or shock) can develop under certain circumstances 
after a finite distance of propagation. Before a carrier 
shock forms, or before its onset is halted by (e.g.) dis- 
persion or the nonlinear response time, the carrier-wave 
undergoes self-steepening. This phenomenon received lit- 
tle attention from the optics community for more than 
30 years, until it was revisited in the 1990s by Flesch et 
al. flip] , [l6| , who performed finite-difference time-domain 
(FDTD) simulations of the process; more recently, the 
role of dispersion was examined in more detail JlJ] . 

The theory of carrier shocking and self-steepening is an 
important way of determining the profiles of pulses with a 
maximized field gradient, and so is needed when consider- 
ing the limits of a gradient-gated HHG scheme. Although 
a fertile field for theoretical work, CSS pulses are difficult 
to realize experimentally because of the deleterious effect 
of dispersion. Nevertheless, the theory provides impor- 
tant information when maximizing the gradient-gating 
effect in HHG. 



A. General case 

We use the method of characteristics (MOC) |l^] to 
predict the shocking distance for a pulse in a dispersion- 
less medium containing an arbitrary combination of non- 
linear terms. The theory is based on the one-dimensional 
(ID), sourceless, plane-polarized Maxwell's equations for 
a field propagating in the z direction, and a material re- 
sponse characterized by the electrical displacement field 



(m) 



E" 
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m> 1 



where X refers to the m-th order nonlinear susceptibil- 
ity of the medium, which is assumed to be instantaneous. 
This extends the model used by Kinsler since it al- 
lows for combinations of all orders of nonlinearity. Prop- 
agating pulses with a field profile E(t) using this model 
leads to self-steepening of the carrier-wave. Eventually, 
a localized infinite field gradient occurs on the wave form 
at a distance 
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Some exotic features of this prediction are discussed in 
JTsj ; a longer derivation is given in appendix |A|. 

To verify the predictions of eqn. (|2j), we used the pseu- 
dospectral spatial domain (PSSD) method jl9| to prop- 
agate pulses for various orders of nonlinearity. The nu- 
merical shocks were detected using the local discontinuity 



detection (LDD) method jl2), which is useful for shock 
detection within a discretized system. We also used a hy- 
brid propagation technique, combinin g d irectional vari- 
ables |2(J and a wideband envelope [El), to check the 
results of the steepened and shocked wave forms. In all 
cases, the LDD (numerical) shocking distances obtained 
from our simulations were in good agreement with those 
predicted by eqn. (|J). 
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1: Self-steepening at the LDD shocking distance for (a) 
X v ~'> (b) X^ 3 \ an< l ( c ) X^ 4 ' nonlinearities. Although a mate- 
rial with only a x' 4 ' nonlinearity is not realistic, its inclusion 
helps illustrate some important features. We can test inver- 
sion symmetry by replacing E with — E. The initial pulse has 
E — Eo s'm(ujt + 0i)sech(O.28wt/r), with a wavelength of 800 

nm, <j>i — 0, and r = 1; no = L The nonlinear strength is 
x (n) E n-i = o Q2 

Fig. |l| illustrates typical CSS pulses for low-order non- 
linearities, and displays some striking properties. In par- 
ticular, figs, [l] (a) and|l| (c) (for and X*- 4 ') both have 
an underlying sawtooth shape, where the unit of repeti- 
tion is a whole cycle. This is because both even and odd 
harmonics are present in the spectrum, and the inversion 
symmetry of the initial pulse is broken; positive lobes 
of E(t) lean to later times, and negative lobes to earlier 
times, since (for x^ 2 "*) the change in characteristic veloc- 
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ity is proportional to dE/dt, and (for x^ 4 "*) to dE 3 /dt. 
In contrast, fig. [j] (b) (for x'" 3 '') shows both positive and 
negative lobes leaning to later times, since the change in 
characteristic velocity is proportional to dE 2 /dt. Because 
odd nonlinearities produce a cascade containing only odd 
harmonics, every half cycle is similar and the pulses pre- 
serve inversion symmetry. Another noteworthy feature is 
that, for higher-order nonlinearities, the steepening be- 
comes increasingly localized to within the most intense 
carrier oscillation(s). 



B. Special case: x' 2 ^ CSS pulses 

We now present a brief review of self-steepening and 
shocking for x^ nonlinearities. We focus on this case 
because it has not previously been addressed and because 
of the relevance to HHG discussed later in Section IV. Of 



course, nonlinearities of higher order than third also exist 
(see, e.g., (2gj), but usually they are much less important. 
The x case has already been extensively addressed in 
the literature fO, g @. For a x' 2 ' nonlinearity 
on its own, eqn. (|) indicates that the shocking distance 
depends primarily on the inverse field gradient so that 
S 2 oc Mm\dE/dt\-\ 

A profile for a (nearly shocked) x^ CSS pulse in a 
dispersionless medium can be seen in fig. ^(a), with its 
characteristic sawtooth profile. It exhibits two distinct 
regions; one containing a ramp- like smooth progression 
from negative to positive field values taking up the ma- 
jority of the cycle, and the second with a rapid positive 
to negative transition over the remaining small fraction 
of the cycle. If we plotted the wave form for a continu- 
ous wave field, rather than a pulse, we would see the field 
evolving towards a sawtooth-like profile (e.g. as seen in 

The build-up of harmonic content during the reshaping 
and steepening of the time domain profiles is vividly il- 
lustrated in fig. |^. Note that the spectra differ from those 
produced under x*- 3 ^ self-steepening where only odd har- 
monics are produced. The nature of the harmonic cas- 
cade in fig. |2| indicates that, if it were possible to achieve 
strong self-steepening experimentally, this would in itself 
be a route to HHG. We will return to this point in the 
next section. 

As for the x^ cas e |jl2), the process is also sen- 
sitive to the pulse length and the carrier envelope phase 
(CEP) offset; theseparameters are also important when 
driving HHG. Fig. |l|(a) shows that for a CEP of fa = 0, 
the carrier has a positive gradient at the center of the 
pulse. However, steepening occurs for negative gradi- 
ents in the x^ case, and hence the high gradients occur 
near the zeros half a cycle away on either side. Since 
these are reduced by envelope fall-off, S2 increases as the 
pulse width is reduced for fa = 0. In contrast, when 
fa = 7r, steepening occurs at the center of the pulse, and 
S2 is then independent of pulse duration. The depen- 
dence of S2 on CEP for different pulse widths is shown 




FIG. 2: (Color online) Logarithmic plot of E(ui) 2 every half 
micron, approaching 5*2 in a material. The spectra flatten 
out as increasing amounts of ever higher-order harmonics are 
generated as the pulse approaches the shock. The pulse and 
material parameters are as for fig. |l[ Each ordinate division 
represents five orders of magnitude. 



in Fig. m This can be compared with the more com- 
plicated CEP dependence seen in fig. 4 of fl2|] , where 
S 3 oc Min \dE 2 /dtr 1 . 
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FIG. 3: MOC shocking distances S2 as a function of CEP (j), 
allowing for different pulse widths. In the long pulse (r = 28) 
case, a shallow minimum centered around 0i — n/2 is just 
visible, caused by the term in the square root of eqn. (Q). 
Apart from the varying pulse length, the pulse and material 
parameters are as for fig. hi. 



III. GENERATING CSS PULSES 

The theory of Section [n] was based on the ideal case of 
a dispersionless medium with instantaneous nonlinearity. 
The question of whether sufficiently self-steepened pulses 
could be generated in real nonlinear materials has been 
addressed for the x^ case in O, [H3|, where considera- 
tions such as linear dispersion [ ]f2| and the strength and 
response time of the nonlinearities were evaluated. Simi- 
lar issues arise in the x case, and so we will not repeat 
them here. We note that for many x^ materials, the 
polarization of the generated harmonics is different from 
that of the driving field, bringing in the further compli- 
cation of birefringence. However, materials do exist that 
support processes where the polarizations are the same 
(e.g. e + e — ► e in lithium niobate [[23]), which could (in 
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principle) generate a harmonic cascade in a single polar- 
ization. 

In practice, the limitations imposed by dispersion and 
material damage mean that generating strongly self- 
steepened pulses would be an extremely challenging task, 
although some success might be achieved if relatively 
weak self-steepening were sufficient. We therefore con- 
sider the possibility of synthesizing CSS-like pulses by 
combining a sub-set of harmonic components using the 
phase offsets and amplitudes obtained from numerically- 
generated CSS profiles. It turns out to be possible to ap- 
proximate an d X*' 3 '' CSS-like wave forms using field 
strengths and phases for the n-th harmonic (n > 1) based 
on the formulae 



(fin = ncf>i + (n — 1) 7r /N. 



(3) 
(4) 



where n(> 1) is the harmonic order, N = 2 or 3 for the 
%( 2 ) and x cases respectively, and F is a fitting factor. 
With F = 4, these equations give a reasonable match to 
the CSS wave form with the spectrum fourth from the 
top in fig. |2| which is close to shocking. We note that 
higher-order harmonics comprise only a small fraction of 
the total CSS pulse energy. However, only a small num- 
ber of harmonics would be used in a synthesized pulse. 
Quite apart from the experimental complexity involved 
in combining a large number of harmonics, the HHG ex- 
periment would be pointless if the high-order harmonics 
were available already in the driving pulse! As we will 
show in Section IV, sufficient gradient enhancement can 
be achieved with a driving pulse containing only two or 
three additional harmonic components. 



IV. HHG AND CSS PULSES 

As described in the Introduction, the bursts of XUV 
radiation produced in the electron ionization-recollision 
process that occurs when an intense few-cycle laser pulse 
interacts with an atom, contain a wide range of high- 
order harmonics. HHG experiments aiming to produce 
isolated attosecond bursts are typically based on few- 
cycle driving pulses, but a number of sophisticated tech- 
niques have also been deployed, including the use of 
two-color driving fields |^, |||lt| , polarization-gating 
§, |, |, §, and chirp control @. 

A typical HHG spectrum from an atomic gas consists 
of a broad plateau of harmonics extending to high or- 
ders, which falls off rapidly above a cut-off energy £. For 
monochromatic driving fields jj], |[ , 



£ = Ip + 3.17U P , 



(5) 



where Ip is the ionization potential of the atom, Up = 
Ef I (2m e u) 2 is the ponderomotive potential of the laser 
field, and E\ is the maximum field strength. The atom 
used for all the simulations in this paper was neon, which 



has a Up of 2I.6eV. Eqn. (|5|) suggests that £ can be in- 
creased either by increasing E\ or by using an atomic 
species with a higher Ip. However, raising E\ also en- 
hances the ionization rate, and this in turn can deplete 
the ground state, causing defocusing of the laser beam 
and a reduction in HHG. 

For few-cycle driving pulses, the parameters naturally 
change from one half-cycle to the next. This feature 
is highlighted in Fig. [|, where a wavelet transform dis- 
plays the variation of the harmonic spectrum for an 8-fs 
driving pulse, calculated using the ID time-dependent 
Schrodinger equation (TDSE) model p6| . For pulsed 
(or multi-color) excitation, the strong field approxima- 
tion (SFA) |27]] predicts the spectral cut-off energy to be 
the maximum value of the function p8| 

E(t u t r ) = I P + ^[eA(t r )+p{U,t r )} 2 , (6) 

where ti is the ionization time, t r the recollision time, 
A(t) is the vector potential, and p is the asymptotic mo- 
mentum of the electron. Eqn. @ reproduces eqn. (||) 
in the case of a monochromatic field. For electrons that 
both ionize and recombine at the position of the core, 



p(U,U) 



tr t{ 



A(t)dt. 



(7) 



Although a general interpretation of eqn. (^|) is com- 
plicated, the electron trajectories that are of most impor- 
tance for HHG are those for which the initial momentum 
of the electron (p(ti, t r ) + eA{ti)) is close to zero. In this 
case, p(ti,t r ) = —eA(ti), so that eqn. (||) simplifies to 
£ = I P + [eA(t r ) - eA(ti)} 2 /2. Usually A[ti) is small, so 
it is broadly true to say that £ depends on the peak value 
of A(t r ). And since the second time derivative of A and, 
concomitantly, the first time derivative of E tend to fol- 
low A, the implication is that maximizing the gradient of 
the field will raise the allowed peak value of A(t r ). Since 
CSS-like pulses have strongly steepened field gradients, 
they should therefore produce HHG with a raised cut-off 
energy £ and an increased number of harmonics. 



A. HHG using synthesized driving pulses 

To create a two-color driving field, a second harmonic 
component is added so that 

E{t) = Eico8(ut + <t>i)+E 2 co$(2ut + <t> 2 ), (8) 

where we will assume that cf>i = for simplicity. The 
addition of the second harmonic term causes inversion 
symmetry of the profile to be lost. If 02 = 0, the pos- 
itive and negative lobes of the field have different peak 
amplitudes, leading to different levels of tunnel ioniza- 
tion on alternate half cycles, and associated variations 
in the intensity and cut-off of the HHG spectrum. On 
the other hand, when <p2 — tt/2, it is the positive and 



4 



Radnor-CKN-2008 



CHI2HHG-II 



Dr.Paul.Kinslcr@physics.org 



100 




Time [fs] 



FIG. 4: (Color online) Wavelet transform of the HHG from 
a normal driving pulse with A =800 nm, a 8-fs FWHM with 
peak intensity 5 x 10 14 W/cm 2 . Since the driving pulse retains 
inversion symmetry, high-order harmonics with a similar dis- 
tribution are generated every half-cycle. Note that the centre 
of theoulse produces the bulk of the HHG, which [following 
eqn. (§)] depends on E 2 . The ID TDSE model of HHG was 
used to generate this figure. 




Time [fs] 



FIG. 5: (Color online) Wavelet transform of the HHG from 
a synthesized x' 2 ' CSS-like driving pulse using amplitudes 
and phases from eqns. (^|) and Q with F — 4 and spectral 
contributions above the sixth harmonic removed. The energy 
and other relevant characteristics of the driving pulse are same 
as for fig. ^. It clearly shows the alternating emission from 
each half cycle of the steepened pulse. The ID TDSE model 
of the HHG process was used to generate this figure. 



negative slopes of E that differ, along with the positive 
and negative excursions of A and d 2 A/dt 2 . Once again, 
the effect is to maximize £ every other half-cycle [ p)| , 
but the effect is now substantially stronger than in the 
02 — case. 

More complex fields can be synthesized by adding fur- 
ther harmonics, with each additional component provid- 
ing a significant boost to the field gradient. However, 
as we shall see in subsection IV C below, the benefits 



are eventually subject to the law of diminishing returns. 
Pulse profiles containing even as well as odd harmonics 
invariably exhibit characteristics that alternate from one 
half-cycle to the next. This feature is clearly evident in 
Fig. H where the wavelet transform of the HHG gener- 
ated by a CSS-like pulse spectrally truncated above 
the sixth harmonic is displayed; details are given in the 
caption. The figure should be compared with Fig. |J. If a 
high-pass filter (such as a molybdenum-silicon multilayer 
mirror) is applied to the signal of Fig. |5|, the lower peaks 
can be removed, leaving the highest peak(s) in greater 
isolation. As such, these filters are useful for the genera- 
tion of a single XUV burst, or when creating burst trains 
for stroboscopic imaging [p9[. 



B. Phasing of the harmonic components 

Increasing the number of harmonics naturally increases 
the experimental complexity, and arranging for their cor- 
rect phasing makes the situation even more complicated. 
It is therefore important to check that CSS-like pulses 
are better for HHG because of their steep field gradi- 
ents, rather than simply because of their greater har- 
monic content. Is the mere presence of many harmonics 
(irrespective of their phase) the crucial thing, or will a 
restricted range of harmonics, correctly phased to maxi- 
mize the gradient, perhaps yield a better result with less 



experimental effort? To address this issue, we compare 
HHG spectra for different types of driving pulse in fig. |^. 
Reading from top to bottom, spectra are shown for 

• a synthesized CSS-like pulse containing third and 
fifth harmonics (labeled x (3) CSS); 

• a synthesized CSS-like pulse containing second, third 

and fourth harmonics (labeled x CSS); 

• a "field-maximized spectrum" (FMS) pulse contain- 

ing identical spectral content to the x^ CSS pulse, 
but with phases <j) n — 0, to create a maximized field 
rather than a maximized gradient; 

• a normal quasi-monochromatic pulse similar to the 

one used to create fig. ||. 

The benefit of a gradient-maximized pulse over both 
the normal and FMS pulses is demonstrated vividly in 
the figure. The highest £ is clearly obtained for the x^ 
CSS pulse. Further, a comparison of x (2) CSS, x (2) FMS, 
and the normal spectra reveals that the CSS pulses have 
the most gradual cut-off. The conclusion is that opti- 
mum performance results from the character of the CSS 
profiles and not merely their harmonic content. 

Fig. shows the temporal profiles of attosecond XUV 
bursts obtained by appropriate filtering of the four HHG 
spectra of fig. ||. In each case, the filtering was adjusted 
to minimize the full width at one quarter maximum, a 
procedure that created the cleanest pulses. It is clear 
that a x CSS driving pulse produces the shortest and 
most isolated XUV bursts, although the intensity is lower 
than in the x^ FMS case where there is one strong and 
one weak burst per cycle. The x CSS pulse produces 
two short XUV bursts per cycle, with a relatively high 
intensity. We note that the oscillations in the tails of the 
X^ CSS spectra in fig.|6|are associated with the presence 
of more than one XUV burst, as seen in Fig. 
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FIG. 6: (Color online) HHG spectra for different driving 
pulses, offset vertically for clarity. Each ordinate division 
represents four orders of magnitude. The pulse energy and 
other relevant characteristics are the same as for fig. kl The 
synthesized y' 2 - 1 FMS and CSS pulses contain second, third, 
and fourth harmonics; the synthesized y 1 - 3 - 1 CSS pulse con- 
tains third and fifth harmonics. Bottom curve to top curve: 
normal (red), y (2) FMS (blue), y (2) CSS (black), and x (3) 
CSS (magenta). 
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FIG. 7: (Color online) XUV burst intensity profiles and 
FWHMs resulting from high-pass-filtered HHG spectra for 
the various CSS driving pulses seen in fig. |(| Bottom curve to 
top curve: normal (red), X (2) FMS (blue), y (2) CSS (black), 
and x (3) CSS (magenta). The curves have been individually 
scaled to aid visibility. 























































































> — e 



300 



240* 



120- 



2 3 4 5 6 7 

Max. harmonic content 

FIG. 8: (Color online) XUV burst durations and intensities 
resulting from high-pass-filtered HHG spectra for synthesized 
X^ 2 ' CSS-like driving pulses, as a function of the maximum 
harmonic component included. The burst durations are de- 
noted with x and use the left-hand scale; the intensities are 
denoted with o and use the right-hand scale. The ID TDSE 
model was used. 



as extra harmonics are added. The FWHM decreases 
roughly linearly up to about the fifth harmonic, but the 
graphs flatten out thereafter. The intensity of the bursts 
drops off rapidly too, although this is not necessarily a 
serious consideration; after all, much intensity is also lost 
in polarization-gating schemes. Different filtering strate- 
gies (e.g. allowing lower frequencies through) can provide 
shorter XUV bursts than shown in fig. || but this usually 
results in stronger satellite pulses. 
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C. Optimum number of harmonics 

As mentioned earlier, moving from standard single- 
color to two-color pumping adds experimental complex- 
ity, and the situation will naturally be exacerbated if 
additional harmonics are included. It is therefore im- 
portant to study how effectively the extra experimental 
investment is reflected in improved HHG performance. 

Figs. H and [| demonstrate what happens when HHG 
is driven by x and X CSS-like pulses as the number 
of harmonic components is increased. In each case, there 
is clearly a strong initial reduction in the burst durations 



FIG. 9: (Color online) XUV burst durations and intensities 
resulting from high-pass-filtered HHG spectra for synthesized 
x' 3 ' CSS-like driving pulses, as a function of the maximum 
harmonic component included. The burst durations are de- 
noted with x and use the left-hand scale; the intensities are 
denoted with o and use the right-hand scale. The ID TDSE 
model was used. 



D. Optimum field profile 

The fact that x CSS pulses succeed in maximizing 
the HHG cut-off energy £ at a given pulse energy is due 
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to the way in which the particular electric field profile 
controls the motion of the ionized electrons. As shown 
in fig. [h], the outward journey and initial deceleration 
of the electrons in group A are determined by the shal- 
low ramp of the field as it rises gently from negative to 
positive. The high return acceleration is managed by the 
high-field region, and the subsequent steep field decrease 
allows the electrons to recollide with the core within mini- 
mal deceleration, whilst simultaneously resetting the field 
amplitude so that the process can repeat. The overall ef- 
fect is to create trajectories that minimize the distance 
traveled, while ensuring the highest possible recollision 
energies along with a high £. The high-energy recol- 
lisions take place largely during the brief high-gradient 
part of the field profile; and give rise to the narrow peaks 
in fig. ||, which then become the short duration XUV 
bursts of fig. after high-pass filtering. 

A genetic algorithm |30|] used in combination with clas- 
sical simulations of electron trajectories, has confirmed 
that the sawtooth-like CSS profiles are uniquely ef- 
ficient for HHG. Using parametrized field profiles, we 
sought to optimize the recollision energy while holding 
both the pulse energy per cycle and the periodicity fixed. 
The algorithm consistently predicted that the optimum 
wave form was a linear ramp starting at a field —E max /2 
and increasing to E rnax . This corresponds to the last 
three quarter-cycles of a sawtooth wave form; indeed in 
fig. we can see that SFA predictions for the CSS- 
enhanced trajectories agree with the classical picture and 
start at about a quarter cycle in. Adding a DC bias or 
sufficiently low-frequency field to a x^ CSS-like profile 
might therefore make fuller use of the whole ramp section 
of the wave form, further increasing £. 

A second family of ionization/ recollision trajectories 
(group B in fig. |l^) are initiated during the high-field 
region of the optical cycle. However, these recollisions 
occur in a region of lower field, produce lower recollision 
energies, and are spread over a longer time interval. They 
lead to more intense but longer HHG emission, with a low 
£, represented by the broader and lower peaks in fig. 

Returning to figs. || and ^, we see that the change 
in burst duration (or intensity) for each additional har- 
monic component is not necessarily straightforward or 
monotonic. This serves to illustrate that that HHG dy- 
namics are complicated, and that the general principle 
of maximizing the field gradient could be optimized by 
further use of a genetic algorithm. Indeed, any HHG 
model could be optimized for a specific outcome - e.g. 
for the shortest-duration XUV burst. It would also allow 
the potential benefits of spectral inputs other than the 
harmonic cascade used here to be evaluated. 



V. CONCLUSION 

We have shown how to optimize HHG using gradient- 
gating of the interaction, and that CSS pulses (or in 
practice their synthesized counterparts) are the most effi- 
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FIG. 10: (Color online) Upper frame: the vector potential 
A{i) corresponding to the field profile shown in the lower 
frame. Note how the peaks in A(t) correspond to the re- 
gions of high field gradient. Lower frame: the trajectories of 
only those high-energy ionized electrons that recollide with 
the nucleus to emit high-order harmonic radiation, calculated 
using the SFA model |27j. The dashed (red) line shows the 
CSS field profile, in this case for the centre of a three-cycle 
pulse with peak intensity 5 x 10 14 W/cm 2 . The CSS-enhanced 
trajectories (starting at A, on the left) recollide over a short 
period of time at high field gradients, giving rise to a brief 
XUV emission with larger cut-off energy. The latter set of 
recolliding trajectories (starting B, right) have a lower energy, 
and recollide over a longer time interval at low field gradients. 



cient way of providing the localized steep field gradients 
needed. An important conclusion is that the gradient- 
gated HHG enhancement can be achieved using as few 
as three or four correctly phased extra harmonic compo- 
nents. The study was assisted by our generalized the- 
ory of optical carrier-wave self-steepening. By describing 
wave forms generated by a nonlinear interaction, the 
theory showed how their advantageous symmetry prop- 
erties could be used to generate shorter, more isolated 
XUV bursts. By using a genetic algorithm with selection 
based on the recollision energy of classical electrons, we 
confirmed that indeed the sawtooth-like CSS pulses 
do have an optimal form that increased the HHG cut-off 
energy £ for any given pulse energy. 

In summary, we have shown that the use of CSS-like 
driving pulses for HHG has the potential to generate 
shorter, more isolated XUV bursts, as well as suggest- 
ing possibilities for other HHG optimization schemes. 



APPENDIX A: MOC THEORY 

Our starting point is the pair of ID sourceless, plane 
polarized Maxwell's equations for a field propagating in 
the z direction, namely: 



dE x 

dz 



dH„ 



dH„ 3D, 



-f^o- 



di 



dz 



Ot 



(Al) 



where E x , H y and D x are the electric, magnetic and elec- 
trical displacement fields respectively. We neglect the 
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tensor nature of the nonlinear coefficients since the po- 
larization of the generated field can be the same as that 
of the incoming field in some situations. Under these 
circumstances, and writing E — E x for convenience, we 
have: 



D = e E + x {1) E- 



(A2) 



m>l 



where x^ m ^ refers to the m-th order nonlinear susceptibil- 
ity of the medium, which is assumed to be instantaneous. 
This leads to a factorizable second order wave equation, 
whose forward propagating part is 



o 



where the characteristic velocity v(E) is, 
v{E) = 



J2 m >l m X {m) E m - 1 



(A3) 



(A4) 



Using the MOC, it is possible to obtain an analytical 
formula for the shocking distance; the technique is used 
for initial value problems in many fields [ fl7| . The objec- 
tive is to turn a first-order partial differential equation 
(PDE) into an ordinary differential equation by chang- 
ing the co-ordinates. In our case, the second-order PDE 



in eqn. ( [A3] ) separates into two first-order PDEs. The 
characteristic curves then follow the field contours, and 
a shock develops where characteristics converge, which 
occurs if a later point travels faster than an earlier point. 
If the characteristics diverge, no shock develops. 

If one characteristic travels slower than another, they 

, 2 

will intersect after some distance i, since ^ ~ T = 2L.. 

at t ' 



Substituting eqn. (A.4) into this, we arrive at the gener- 
alized formula for the distance at which characteristics 
converge: 



2c k + J2QX (q) E q - 1 

v ?>i 



E 

m>l 



niX 



(m) 



dE n 



(A5) 



dt 



The shocking distance (S) is now the minimum value of 
L, i.e. S — Min|L|. This minimum can be easily found 
numerically. 
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